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Abstract: The Moderate Resolution Imaging Spectroradiometer (MODIS) provides daily global coverage, but the 10 km 


resolution of its aerosol optical depth (AOD) product is not suitable for studying spatial variability of aerosols in 
urban areas. Recently, a new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm was 
developed for MODIS which provides AOD at 1 km resolution. Using MAIAC data, the relationship between 
MAIAC AOD and PM 2 .5 as measured by the 27 EPA ground monitoring stations was investigated. These results 
were also compared to conventional MODIS 1 0 km AOD retrievals (MOD04) for the same days and locations. The 
coefficients of determination for MOD04 and for MAIAC are R 2 =0.45 and 0.50 respectively, suggested that AOD 
is a reasonably good proxy for PM 2 . 5 ground concentrations. Finally, we studied the relationship between PM 2 5 
and AOD at the intra-urban scale «10 km) in Boston. The fine resolution results indicated spatial variability in 
particle concentration at a sub-1 0 kilometer scale. A local analysis for the Boston area showed that the AOD-PM 2 5 
relationship does not depend on relative humidity and air temperatures below ~7 °C. The correlation improves 
for temperatures above 7 - 16 °C. We found no dependence on the boundary layer height except when the 
former was in the range 250-500 m. Finally, we apply a mixed effects model approach to MAIAC aerosol optical 
depth (AOD) retrievals from MODIS to predict PM 2 5 concentrations within the greater Boston area. With this 
approach we can control for the inherent day-to-day variability in the AOD-PM 2 . 5 relationship, which depends 
on time-varying parameters such as particle optical properties, vertical and diurnal concentration profiles and 
ground surface reflectance. Our results show that the model-predicted PM 2 5 mass concentrations are highly 
correlated with the actual observations (out-of-sample R 2 of 0.86). Therefore, adjustment for the daily variability 
in the AOD-PM 2 . 5 relationship provides a means for obtaining spatially-resolved PM 2 5 concentrations. 
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1. Introduction 

The adverse health effects from particulate matter (PM) 
pollution with aerodynamic diameter ^ 2.5 pm (PM 2 . 5 ) 
must be considered In developing policies to Improve 
air quality [1]. Substantial epidemiologic literature has 
demonstrated that exposure to ambient particulate mat- 
ter (PM) Is associated with Increased morbidity and mor- 
tality, particularly associated with cardiopulmonary dis- 
ease [2, 3] and lung cancer [3, 4]. Since PM Is created 
by various anthropogenic and natural sources with vastly 
different atmospheric residence times, It has a high spatio- 
temporal variability. An accurate assessment of this vari- 
ability Is Important as It leads not only to stronger associ- 
ations between exposure and health but also to deeper un- 
derstanding of epidemiological time-series studies [5, 6 ]. 

Routine measurements of ground-level PM2.5 concentra- 
tions by air quality monitoring networks are of great Im- 
portance In assessing exposures, but their spatial coverage 
has been limited. However, recently It has become clear 
that satellite remote sensing can be an Important tool to 
complement the ground level measurements. The relevant 
satellite-derived parameter Is the aerosol optical depth 
(AOD) which quantifies the extinction of solar radiation 
at a given wavelength due to presence of aerosols In an 
atmospheric column. Because the satellite-derived AOD 
Is a measure of how much light Is absorbed/scattered by 
particles In the column that are affected by ambient con- 
ditions (e.g., variable humidity and consequently variable 
amounts of water adsorbed on particles), while PM2.5 mass 
Is a measure of dry particles near the surface, these two 
parameters are not expected to be strictly correlated. Fur- 
thermore, to be used for air quality applications, Includ- 
ing health studies, the satellite retrieved AOD data (e.g. 
a total column optical measurement) must be converted 
to estimates of PM2.5 concentrations (e.g. a surface-level 
particulate mass measurement). This type of analysis re- 
quires PM2.5-AOD collocated pairs which Itself Is a re- 
strictive requirement [7]. 

Until recently, the main source of global satellite aerosol 
data was the MODIS satellite MOD04 algorithm, which 
provides data at a 10 km resolution. Recently, a new 
Multi-Angle Implementation of Atmospheric Correction 
(MAIAC) algorithm was developed for MODIS which pro- 
vides aerosol Information at 1 km resolution [ 8 , 9]. We 
started with a direct comparison between MOD04 and 
MAIAC retrievals. Toward this end, we conducted a multi- 
year analysis to study the relation of same-day/same lo- 
cation AOD vs PM 2 . 5(2002 - 2008) In the southern part 
of New England. In addition, we conducted a multi-year 
analysis by breaking down AOD vs PM 2.5 regressions by 
season (spring, summer, fall, winter) and by site location. 


Furthermore, we studied the Intra-urban (at scales less 
than 10 km) variability of the relationship between PM2.5 
and AOD for Boston. Finally, we explored whether It was 
possible to obtain accurate estimates of PM2.5 concen- 
trations using a MAIAC AOD retrieval and mixed effects 
model approach (dally adjustment for AOD vs PM2.5 re- 
lationship). Our goal was to show how variability In the 
AOD vs PM2.5 relationship can be captured by a statis- 
tical model during one year of data (January 1 through 
December 31 , 2003). 

1.1. Ground level PM 2 ^observations 

Twenty-four hour-integrated PM 2.5 concentrations were 
measured at 26 U.S. Environmental Protection Agency 
(EPA) PM 2.5 monitoring sites during 2002 - 2008 (Fig- 
ure 1, highlighted by dots and Table 1). These Include 
15 sites from Massachusetts (MA) and 11 sites from Con- 
necticut (CT). Sampling frequency differed by site and In- 
cluded samples collected every day, every third day, and 
every sixth day. Additionally, we used 24 hour-integrated 
PM 2.5 concentrations from the Harvard School of Public 
Health (HSPH) supersite located near downtown Boston. 
Data from this site have been used In a large number of 
epidemiological studies to assess the temporal variability 
of Individual and population exposures In the region. 
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Figure 1 . Study area and EPA monitoring sites for New England 
used for comparison between MOD04 and MAIAC data. 
Area highlighted by box is the Boston urban domain. 
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Table 1 . EPA ground monitoring sites used in this study over New 
England. Boston sites are highlighted in italic font. 


Site ID 

City 

Latitude Longitude 

09-001-0010 

Bridgeport, CT 

41.17 

-73.19 

09-001-0113 

Bridgeport, CT 

41.18 

-73.19 

09-001-1123 

Danbury, CT 

41.40 

-73.44 

09-001-2124 

Stamford, CT 

41.06 

-73.53 

09-001-3005 

Norwalk, CT 

41.11 

-73.41 

09-001-9003 

Westport, CT 

41.12 

-73.34 

09-003-1003 

E. Hartford, CT 

41.78 

-72.63 

09-003-1018 

Hartford, CT 

41.76 

-72.67 

09-009-0018 

New Haven, CT 

41.29 

-72.90 

09-009-0026 

New Haven, CT 

41.29 

-72.89 

09-009-1123 

New Haven, CT 

41.31 

-72.92 

09-009-2008 

New Haven, CT 

41.33 

-72.92 

09-009-2123 

Waterbury, CT 

41.55 

-73.04 

09-009-8003 

W. Haven, CT 

41.28 

-72.96 

09-011-3002 

Norwich, CT 

41.52 

-72.08 

25-005-1004 

Fall River, MA 

41.68 

-71.17 

25-009-2006 

Lynn, MA 

42.47 

-70.97 

25-009-5005 

Haverhill, MA 

42.77 

-71.10 

25-013-0008 

Chicopee, MA 

42.19 

-72.56 

25-013-0016 

Springfield, MA 

42.11 

-72.59 

25-013-2009 

Springfield, MA 

42.11 

-72.60 

25-023-0004 

Brockton, MA 

42.08 

-71.01 

25-025-0027 

Boston, MA 

42.37 

-71.06 

25-025-0042 

Boston, MA 

42.33 

-71.08 

25-025-0043 

Boston, MA 

42.36 

-71.05 

25-025-0002 

Boston, MA 

42.35 

-71.10 

25-027-0020 

Worcester, MA 

42.27 

-71.80 

Harvard supersite 

Boston, MA 

42.34 

-71.10 


1.2. Satellite data 

A new algorithm MAIAC [8] has been developed to process 
MODIS data. MAIAC retrieves aerosol parameters over 
land at 1 km resolution simultaneously with parameters 
of a surface bidirectional reflectance distribution function 
(BRDF). This Is accomplished by using the time series 
of MODIS measurements and simultaneous processing of 
a group of pixels. The MAIAC algorithm ensures that 
the number of measurements exceeds the number of un- 
knowns, a necessary condition for solving an Inverse prob- 
lem that does not regulre the assumptions typically used 
by current operational algorithms. The MODIS time se- 
ries accumulation also provides multi-angle coverage for 
every surface grid cell, which Is regulred for the BRDF 
retrievals from MODIS data. The aerosol parameters In- 
clude optical depth (total aerosol) and fine mode frac- 
tion. Following the MODIS operational aerosol algorithm 


(MOD04) [9], models for the fine and coarse aerosol frac- 
tions are specified regionally based on the climatology 
of the Aerosol RObotic NETwork (AERONET) [10] sun- 
photometer data. AERONET validation over the conti- 
nental USA showed that the MAIAC and MOD04 algo- 
rithms have a similar accuracy over dark and vegetated 
surfaces, but also showed that MAIAC generally Improves 
accuracy over brighter surfaces, Including most urban ar- 
eas [8]. The Improved accuracy of MAIAC resulted from 
the explicit surface characterization method, In contrast 
to the empirical surface parameterization approach, uti- 
lized In the MOD04 algorithm. Furthermore, MAIAC In- 
corporates a cloud mask (CM) algorithm based on spatio- 
temporal analysis which augments traditional pixel-level 
cloud detection technlgues [11]. In addition, the residual 
contamination by clouds and cloud shadows was reduced 
by discarding 2 pixels adjacent to detected clouds. 

In addition to MAIAC data we used dally MODIS Level 
2 (MOD04) Collection 5.1 Aerosol data from the Terra 
platform that are produced at a spatial resolution of 
10x10 km 2 (at nadir). More details about the MODIS 
AOD retrieval are reported In [8, 12]. We conducted a 
comparative analysis of AOD between MAIAC and the re- 
spective MOD04 product. It Is Important to mention that 
the MOD04 product Is reported for an area of 20 by 20, 
500 m pixels In the swath format. This area corresponds 
to spatial resolution of 10x10 km 2 at nadir, however It 
grows with the scan angle reaching -20x40 km 2 at the 
edge of scan due to the respective growth of the MODIS 
pixel footprint by a factor of -2x4. Conversely, MAIAC 
provides a uniform 1 km grldded resolution at selected 
projection regardless of the scan angle. This means that 
the MAIAC product Is under-sampled by a factor of 4 at 
nadir, considering maximal available spatial Information 
from 500 m pixels, and Is oversampled by a factor of 2 
at the edge of scan. In this regard, MOD04 data are 
always under-sampled by a factor of 400. In order to 
perform a direct MOD04-MAIAC comparison, the area of 
each MOD04 pixel was approximated by a polygon, and 
all MAIAC 1 km data fitting this area were averaged. 

1.3. Meteorological data 

All meteorological variables used In the analysis (tem- 
perature, boundary layer and relative humidity) were ob- 
tained through the National Climatic Data Center (NCDC) 
(NCDC, 2010). Only continuous operating stations with 
dally data running from 2000 to 2008 were used. In ad- 
dition, we used meteorological data from Boston Logan 
airport. Grid cells were matched to the closest weather 
station for meteorological variables. 


High resolution aerosol data from MODIS satellite for urban air quality studies 


1.4. Data analyses 

We investigated the associations between AOD and PM 2.5 
daily measurements at the sampling sites for the years 
2002-2008. We first made a direct comparison between 
MOD04 and MAIAC retrievals, with a multi-year analy- 
sis of AOD vs PM 2.5 for the same days (2002-2008) and 
locations (27 EPA monitoring stations) in New England. 
Using the same data we performed AOD vs PM 2.5 regres- 
sion analysis by season (spring, summer, fall, and winter) 
for each of the three regions. In addition, we conducted 
AOD vs PM 2.5 regression analysis by site location. 

Next, we analyzed the intra-urban variability in AOD vs 
PM 2.5 relationship inside a 10 x 10 km 2 area of greater 
Boston containing four EPA ground monitors and the Har- 
vard Supersite. To take into account the variability of 
the sampling frequency of the EPA stations, only days 
with at least three available AOD-PM 2.5 pairs inside a 
10x10 km 2 box were selected (Figure 1, highlighted by 
the polygon). There were a total of 304 days with 3 - 
5 observations. In addition, we analyzed PM 2.5 vs AOD 
association in relation to PM 2.5 particle levels, temper- 
ature, boundary layer height and relative humidity. Fi- 
nally, we explored whether it is possible to obtain ac- 
curate estimates of PM 2.5 concentrations using MAIAC 
AOD retrievals and statistical modeling with a resolution 
of 1 xl km 2 conducted for Boston, MA in the northeastern 
part of U.S. 

1 .5. Mixed effects model approach 

In this study we used a mixed-effects model approach that 
accounts for day-to-day variability. A basic assumption is 
that the relationship varies daily because it depends on 
time-varying parameters such as relative humidity, PM 2.5 
vertical and diurnal concentration profiles, PM 2.5 optical 
properties and surface reflectance. In a recent paper, we 
showed that the mixed effects model approach provides 
higher accuracy and precision in predicting PM 2.5 con- 
centrations based on the MODIS AOD dataset than a 
simple regression model [13, 14]. In the present study we 
use this model approach to predict PM 2.5 concentrations 
based on MAIAC AOD retrievals. Conseyuently, guan- 
titative relationships between PM 2.5 concentrations mea- 
sured at the 27 PM 2.5 monitoring sites and AOD values 
in their corresponding grid cells were derived. We used 
the following mixed effects model with random intercepts 
and slopes (Eg. 1): 

PMij = (a+Uj)+[(p'\+Vj)xAOD i j]+Si+£ij(UjVj) ~ [(oo),Tp] 

( 1 ) 

where PM ty is the PM2.5 concentration at a spatial site 
i on day j; AOD /y is the AOD value in the grid cell cor- 


responding to site i on day j; a and u y are the fixed and 
random intercepts, respectively; fa and v y are the fixed and 
random slopes, respectively; Si ~N(0, as2) is the random 
intercept of site i; cij ~N(0, a2) is the error term at site 
i on a day j; and L/3 is the variance-covariance matrix 
for the random effects. The AOD fixed effect in the model 
(Eg. 1 ) accounts for the effect of AOD on PM 2 . 5 , which was 
the same for all study days. The AOD random effects ex- 
plain the daily variability in the PM 2 . 5 -AOD relationship. 
The solution of the mixed model equations is a maximum 
likelihood, a form of estimation that accounts for the pa- 
rameters in the fixed-effects structure of the model to re- 
duce the bias in the covariance parameter estimates [15]. 
Currently, this is the method implemented for the SAS 
statistical software package (proc mixed). 

Finally, PM2.5 concentrations for each grid cell on a day 
j were estimated using the corresponding AOD values as 
follows: 

PMij = (a + uj) + [(fa + Vj) x AODij] + c, y (2) 

The fixed and random intercepts, and the fixed and random 
slopes for each study day were derived previously from 
Eg. 1. Note that the random estimates for the site term 
were excluded. In this way AOD values were unbiased 
and representative of their corresponding grid cell. 

We use a cross-validation (CV) approach to evaluate the 
ability of the model to predict PM2.5 concentrations for 
each pixel in the study area. Thus, the dataset was re- 
peatedly randomly divided into 90% (calibration) and 10^6 
(held-out test) splits. We applied the fitted calibration 
model to estimate PM2.5 for the held-out test set. This 
"out-of-sample" process was repeated ten times to calcu- 
late the cross-validated (CV) R 2 values. 

2. Results 

2.1 . Direct comparison with MOD04 retrieval 

This section studies the subset of MOD04/MAIAC data 
for days when both products are available for a given EPA 
site. Figure 2 shows the direct comparison between PM 2.5 
and AOD for MOD04 (10 km) and MAIAC (1 km) for the 
same days and locations (2002-2008) in New England (27 
locations, N=2310, p<0.0001). The coefficients of deter- 
mination (R 2 ) for MOD04 and MAIAC are 0.45 and 0.50 
respectively, suggesting that AOD is a reasonably good 
proxy for PM 2.5 ground concentrations. In other words, 
near-surface PM 2.5 concentrations do not reflect the total 
AOD column values. 

The previous research has shown that the PM2.5 vs AOD 
relationship varies seasonally and by location [16]. Ta- 
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ble 2 presents a multi-year, seasonal (spring, summer, 
fall, winter) comparison between MOD04 and MAIAC. Al- 
though MAIAC shows intercepts that are lower than those 
for MOD04, for 8 year of measurements, slopes for both 
retrievals are similar, with the range of slopes between 7 - 
8 pg/m 3 in winter and 26 - 31 pg/m 3 in spring, summer 
and fall. The slight improvement in correlation is related 
to the finer resolution of MAIAC with its better correspon- 
dence between the monitoring site and the respective grid 
cell size, and better performance over bright urban areas. 
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Figure 3. Frequency distribution of daily AOD vs PM 2.5 correlations. 
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Figure 2. Comparison between PM 2.5 and AOD for MODIS 10 km 
(MOD04) and MAIAC 1 km for the same days and 
locations (2002-2008) in New England (27 locations, 
N=2310). The solid line represents the regression line, 
and the dashed line displays the 1 :1 line. 


Table 2. Seasonal comparison between coarse MOD04 AOD 1 0 km 
and fine resolution MAIAC 1 km AOD for the same days and 
locations. 


Data Source 

Statistics 

Summer 

Fall 

Winter Spring 


N 

786 

886 

74 

564 


Intercept 

8.15 

7.26 

7.41 

5.38 

MOD04 

Slope 

26.3 

28.9 

7.7 

26.2 


p-value 

<.0001 

<.0001 

0.285 

<.0001 


R 2 

0.45 

0.30 

0.002 

0.35 


N 

786 

886 

74 

564 


Intercept 

6.08 

5.8 

7.04 

3.56 

MAIAC 

Slope 

25.3 

28.9 

8.04 

31.9 


p-value 

<.0001 

<.0001 

0.226 

<.0001 


R 2 

0.50 

0.35 

0.007 

0.41 


Figure 3 shows the frequency distribution of the correla- 
tion coefficient between PM 2 .5 and AOD by date for 2002 - 
2008, for the same days and sites. In general, both re- 
trievals provide a similar accuracy. Importantly, as can be 
also seen, the AOD vs PM 2.5 relationship changes by date 
for both, MOD04 and MAIAC indicating a clear temporal 
variation in the association between both parameters. 
While Figure 3 shows the daily variability in AOD vs 


PM 2.5 relationship, Figure 4 presents correlation coeffi- 
cients between PM 2.5 and AOD by the EPA site location 
for 2002-2008 (the same days were used for MAIAC and 
MOD04, p<0.005 for all sites). In general, both retrievals 
provide similar results. Note that the range of correla- 
tions for both retrievals across the sites is substantial, 
which most likely reflects the local meteorological condi- 
tions and spatial homogeneity of PM 2 . 5 , namely how well 
the local PM 2.5 measurement can be generalized to the 
larger footprint of the AOD pixel. 
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Figure 4. Correlation coefficient between PM 2 .sand AOD by EPA 
site location for 2002-2008: dashed line indicates corre- 
lation for MAIAC and solid line for MOD04. 


2.2. Intra-urban variability in the AOD vs 
PM 2 5 relationship 

The high resolution AOD potentially carries information 
about local-scale variability, which is especially impor- 
tant for highly populated urban areas. We define the lo- 
cal variability as the variability in daily averaged PM2.5 
values among different EPA stations (from 3 to 5) in the 
10x10 km 2 box. Figure 5 shows an intra-urban AOD- 
PM2.5 correlation in Boston with R 2 =0.38 for the entire 
study period 2002 - 2008 (N=304 days). Each PM2.5 
value for a given station in Figure 5 represents a daily 
averaged value. 
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AOD 


Figure 5. Intra-urban scatter plot between PM 2 .5 at EPA monitoring 
stations (with at least 3 ground PM 2.5 measurements) and 
1 km MAIAC AOD measured over the entire period 2002- 
2008 for Boston (N=304 days). The solid line represents 
the regression line, and the dashed line displays the 1:1 
line. 


Daily averaged concentrations ofPM 2 5 




Figure 6. Correlation coefficient as a function of spatial variability in 
PM2.5 concentrations. 


temperatures with shallow BLH are associated with win- 
ter conditions when average PM 2.5 and AOD are low and 
one cannot expect good AOD-PM 2.5 correlation. 


First, we explored an intra-urban AOD-PM 2.5 correlation 
during different pollution days, based on EPA observa- 
tions. Figure 6 shows correlation coefficient by date for 
different levels of daily averaged PM 2.5 concentrations. 
As expected, at low PM 2.5 levels (<5 pg/m 3 ), the distri- 
bution of correlation coefficient is almost uniform in the 
range [—1; 1] indicating low sensitivity of AOD and high 
relative errors of both AOD and PM 2 . 5 . The correlation 
improves at higher PM 2.5 levels, notably for ranges 5-10 
and 10-20 pg/m3. The correlation seems to decrease at 
higher PM levels (>20 pg/m3) which may be an artifact 
of low sampling statistics. 

The intra-urban AOD-PM 2.5 correlation (Figure 5) in- 
cludes the temporal meteorological variability for 2002 - 
2008. It has been shown that under conditions of a well- 
mixed boundary layer with low ambient relative humid- 
ity (RH), the relationship between PM 2.5 and AOD may 
be very robust [17]. With this in mind, in Figure 7a-c 
we studied the influence of temperature, relative humidity 
and boundary layer height (BLH) on daily PM 2.5 vs AOD 
correlations. Our results do not show any dependence of 
daily AOD-PM 2.5 correlations in the Boston metropolitan 
area on RH and air temperature below 7°C. The correla- 
tion improves at higher temperatures, in particular above 
7°C typical of late spring-summer-early fall conditions. 
Figure 7b also shows an improved correlation for the mod- 
erate boundary layer heights of 250 - 500 m, and a poorer 
correlation outside of this range. Some of these results 
have a clear physical explanation: for example, the low 


2.3. Prediction of PM 2 5 concentrations using 
mixed effect model approach 

Since these relationships between AOD and PM2.5 mea- 
surements vary daily, mixed effects models were used to 
allow for the regression intercepts and slopes to vary 
daily. The fixed effects of the AOD intercept and slope 
were statistically significant: a =9.3 (p<0.0001) and 
ft =17.2 (p<0.0001) respectively. The fixed effects of 
spatial and temporal predictors were also significant. In 
addition, the random slopes for AOD by day, and by day 
and region were both significant (p<0.0001). Figure 8 
shows the daily variation of random AOD intercepts and 
slopes. Note that these results support the findings that 
because the parameters influencing the relationship be- 
tween PM2.5 and AOD vary from day to day within a given 
domain, it is necessary to adjust for this daily variability. 

The measured and predicted PM 2.5 concentrations in the 
cross validation (CV) model are shown in Figure 9. As can 
be seen, the CV mixed effects model performed quite well. 
The CV test resulted in a R 2 value of 0.90 and slope of 
0.91, indicating a good agreement between the measured 
and predicted concentrations. 

Finally, the model was applied for retrieval days to ex- 
plore the pattern of PM2.5 concentrations on a daily basis. 
June 25th was selected as the high pollution event to an- 
alyze the predicted PM2.5 concentrations resulting from 
the mixed effect model approach. During this day, the av- 
erage concentrations on EPA monitoring stations in the 
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Temperature Boundary layer height (m) 



Figure 7. Correlation coefficient by date conditioned at several meteorological parameters: a) temperature, b) boundary layer height, and c) 
relative humidity. 



Slope 



Intercept 


Figure 8. Frequency distribution of the random intercepts and 
slopes. 


New England area ranged between 25-42 pg/m3. During 
20-24th of June 2003, a forest fire in the Quebec province, 
Canada, brought smoke pollution into Massachusetts, NE. 
Figure 10 presents a MODIS Level IB true-color image 



Figure 9. Measured vs predicted PM 2 .5 concentrations for the cross- 
validation model (left) and the test set (right). The solid 
line represents the regression line, and the dashed line 
displays the 1 : 1 line. 


and shows long range transport of thick haze from fires 
in the Quebec province, Canada on June 24, 2003 (A) 
and June 25, 2003 (B) (left, red arrow) along the east- 
ern seaboard of the U.S. Figure 9 C and D highlight high 
AOD values at 550 nm (right) which are elongated with 
the plum. 

In Figure 11 we show the spatial pattern of PM 2 .5 con- 
centrations resulting from the mixed effect model on June 
25, 2003. As can be seen, the spatial concentration pat- 
terns are different across the domain and are highest in 
Boston due to high pollution transport from forest fires in 
Canada. Importantlg, as shown in Figure 11, the variabil- 
itg in PM 2 .5 concentrations at fine scale resolution can be 
high even during high pollution events when the contrast 
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Figure 1 0. Maps A and B: MODIS Level 1 B true-color image shows 
long range transport of thick haze from fires in the Que- 
bec province, Canada on June 24, 2003 (A) and June 
25, 2003 (B) (left, red arrow) along the eastern seaboard 
of the U.S. Maps C and D: High AOD values at 550 nm 
(right) are elongated with the plume. 



Figure 11. PM 2 .5 concentrations modeled by mixed effects ap- 
proach for June 25, 2003. 


in pollution levels between geographically adjacent areas 
is not expected to be high. Compare dashed line 1 on 
Figure 11 : for the area of -20 km, variability in parti- 
cle concentration ~4 pg/m 3 that can be captured by the 
model. Furthermore, the fine resolution results indicated 
spatial variability at a sub-10 kilometer scale. 

3. Conclusions 

This paper analyzed how the spatial resolution of the 
AOD product affected the correlation between satellite- 
retrieved AOD and ground based PM2.5 concentrations 
using 7 years of MODIS Terra observations over the 
southern part of New England. A direct comparison be- 
tween coarse MOD 04 10 km AOD and high resolution 
MAIAC 1 km AOD for all collocated AOD-PM2.5 pairs for 
the same days and locations showed that although both 
retrievals provide reasonable results, MAIAC was found to 
provide a slightly better correlation. Furthermore, there is 
clear temporal variation in the association between AOD 
and PM2.5. Importantly, a local analysis for Boston area 
showed that AOD-PM2.5 relationship does not depend on 
RH and air temperatures below ~ 7 °C. The correlation im- 
proves for temperatures above 7 - 1 6 °C. We found a poorer 
correlation between AOD and PM2.5 on days with very low 
or very high boundary layer height. 

From the epidemiological and exposure assessment point 
of view, it is of high importance to have information about 
the spatial variability of the exposures in the city. Sev- 
eral studies published in the last 3 years have shown that 
high spatial resolution is essential to detect spatial vari- 
ability in PM levels [ 18 ] and in aerosol loadings at re- 
gional and at a sub-10 km scale (e.g. intra-urban do- 
main) [ 19 , 20 ]. Our study using MAIAC data and mixed 
effect approach showed high accuracy in the New Eng- 
land domain thereby indicating that our model based on 
MAIAC data can be used to investigate the intra-urban 
exposure contrasts in PM2.5 levels. 

Our results show a daily adjustment using a mixed ef- 
fects model approach effectively controls the combined 
effects of many parameters that can influence the daily 
variability in the AOD-PM2.5 relationship. This implies 
that within a given region, the types of aerosols may be 
more homogeneous and the height of the boundary layer 
and humidity may be more uniform, making the relation- 
ship between AOD to PM2.5 less variable. Therefore, the 
proposed method has the advantage that it can easily be 
applied to other regions by taking into account the con- 
ditions prevailing in each region, and adjusting for daily 
variability in the AOD vs PM2.5 relationship. 

Despite promising results, more data need to be processed 
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and analyzed to understand the full potential and lim- 
itations of the high resolution MAIAC AOD product for 
Improving the accuracy In PM2.5 estimations. Next, addi- 
tional parameters can be considered In the model thereby 
Improving the accuracy of PM2.5 estimates. Furthermore, 
In order to further Investigate the strengths and limita- 
tions of using high resolution satellite AOD data for the 
routine modeling of PM2.5 concentrations during high and 
low pollution days we are planning to conduct a compre- 
hensive multi-year study based on the full set of MODIS 
measurements. Next, further Improvements to the MAIAC 
AOD retrieval algorithm would Improve accuracy In PM2.5 
estimation. It Is Important to emphasize though, that the 
Information content of AOD data alone Is limited, and the 
best results may be achieved by combining different data 
sources Including, for example, the aerosol vertical profile 
Information from satellite or ground-based Ildars. 
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